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ABSTRACT 


The Apollo command module aft heat shield, a 
composite laminated torispherical shell, was ideal- 
ized as a network of curved beams. A stiffness 
matrix method was developed in torispherical coor- 
dinates to analyze the structural model for thermal 
and mechanical loadings. Correlation of the analyti- 
cal method was demonstrated with a closed-form 
solution of the Apollo aft heat shield under a cold- 
soak thermal condition. Analyses of the structure 
were also correlated with data from a static test. 
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STIFFNESS MATRIX SOLUTION FOR SHELLS OF REVOLUTION 

POSSESSING VARIABLE THICKNESS HEAT PROTECTION 

By Frederick Jean Stebbins 
Manned Spacecraft Center 


SUMMARY 


A composite laminated torispherical shell was idealized as a network of curved 
beams. A stiffness matrix method was developed to analyze the structural model. 

The method was developed in torispherical coordinates which are amenable to the cur- 
rent family of manned spacecraft. The laminated shell was composed of an ablator, 
a bond, two stainless -steel face sheets, and a 2-inch-thick stainless -steel core. 

A structural model of the Apollo aft heat- shield substructure was analyzed for a 
variety of thermal and mechanical loadings. A correlation was made of the analysis 
with test data and two independent digital computer analyses. An analysis of a sym- 
metrical model of the Apollo aft heat shield under cold- soak conditions was compared 
with a closed-form solution to establish a link with " exact" theory. Selected stresses 
and displacements were used in the comparisons to demonstrate the accuracy of the 
procedure. 


INTRODUCTION 


The high strength-to-weight ratio of shells has been effectively utilized in space- 
craft design. Elements of the common shells of revolution, that is, spheres, cylin- 
ders, tori, and cones, are assembled into spacecraft structural shapes. The Apollo 
spacecraft, for example, is a contiguous assembly of two spherical sections, a trun- 
cated cone and a toroidal section. Hyperthermal conditions during reentry require 
a heat-protective system surrounding the structural shell. For heat protection, the 
United States spacecraft have commonly had a variable thickness ablative material 
which forms a nonstructural char during reentry. The mechanical properties of the 
ablator in the virgin state are highly temperature-dependent. Contemporary space- 
craft design features a stainless-steel sandwich substrate. The ablator and substrate 
form a complex laminated shell which is rotationally and meridionally asymmetric. 

The structural design of spacecraft presents many interrelated problem areas. 
At the time of maximum dynamic loading during the launch phase, the solution to the 
problem involves more than static-load considerations. The dynamic mode shapes of 
the spacecraft must be established and the possibility of buffeting investigated. The 
translunar and reentry flight phases present two opposite extremes in thermal loading. 



The operation of thermally deformed hatches is also a problem. These hatch doors 
must be used during the mission for extravehicular activity by the astronauts and for 
operation of equipment such as the astrosextant. The fiery reentry phase places ex- 
treme requirements on the spacecraft materials; however, the aetual landing creates 
more structural problems than any other phase of the spaceflight mission. These 
structural problems of landing begin with the jolt of the arresting parachutes. Even 
though landing represents the terminal phase of the mission, it is most important that 
the watertight integrity of the spacecraft be maintained since the spacecraft must float 
until it is recovered. 

The requirements imposed by aerodynamic considerations have resulted in geo- 
metric shapes which are not structurally desirable. In the cases of the Mercury, 
Gemini, and Apollo vehicles, the structure of the spacecraft has involved a prominent 
torispherical section. From the viewpoint of the structural engineer, it is desirable 
that the radius of the spherical area be small in order to minimize the water impact 
during landing. However, aerodynamic stability requirements necessitate a relatively 
large spherical radius. The necessary compromise between aerodynamic and struc- 
tural objectives emphasized the need for accurate structural analysis. 

Most of the analytical work in this 
area has been confined to the studies of 
shells of revolution (refs. 1 to 5). The 
objective of the research presented in this 
report was to develop a method of static 
structural analysis for a class of asym- 
metric shells subjected to arbitrary loads 
and temperatures. The thermal loads 
were limited to inplane forces and moments 
having vectors that are tangent to the shell 
midsurface. The thermal gradient may be 
arbitrary in the thickness direction of the 
shell and the surface coordinates. The 
shell may have variable meridional and 
circumferential stiffness properties. Since 
spacecraft, as well as airplanes, have a 
plane of symmetry relative to the flight 
path (fig. 1), the attendant reduction of 
unknowns was made. The application of 
the present method of static structural 
analysis of asymmetric shells under arbi- 
trary loading was made to the Apollo aft 
heat- shield substructure. Theoretical and 
and experimental correlation is demon- 
strated. 
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SYMBOLS 


A AB’ A AB’ a ba^ 
A BA’ A B’B’ a oa| 
B 

h e ,h cp 


E 

f 


BB 


G 

I 


one end of the structural element AB 

transformation matrices of the typical structural element, that is, 
A„ is the transform from end j to end i (ref. 6) 

one end of the structural element AB 
beam width normal to Q or cp direction, in. 

constants to C^g as used in table A-I 

distance from the centroid of composite structure to ith laminate, in. 
modulus of elasticity, psi 

flexibility matrix of the typical structural element for end A fixed 
and end B free 

modulus of rigidity, psi 

4 

moment of inertia, in 

4 

geometrical torsion constant, in 


k 

k AA’ k BB 
k AB’ k BA 
k AB 

P 1 ,P 2’ P 3 
P 4’ P 5’ P 6 


joint stiffness matrix 

stiffness submatrices of the typical structural element for end B 
fixed and end A free, and for end A fixed and end B free, re- 
spectively (ref. 6) 

stiffness submatrices of the coupling actions at end A due to dis- 
placements at end B, or for coupling actions at end B due to 
displacements at end A (ref. 6) 

total stiffness matrix of the member AB 

force in R, cp, and 6 directions, respectively, lb 

moment about R, cp, and 6 axes, respectively, in-lb 
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P A’ P B 
P BA’ P BB 
P A0’ P B0 


R, Rj, Rg 


R 


r 

S 

T 

T 

( 

t 


u l’ u 2’ u 3 


u 4’ u 5’ u 6 


U A’ U B 


U B0 ,U A0 

u. 

X 


a 

y Y y 2 

e 


column vectors of internal actions at end A and at end B respec- 
tively, of a structural element AB 

column vectors of actions at end B due to loads at end A and at 
end B 

column vectors of applied actions at end A and at end B of a 
structural element AB 

group of action components p^ to p g (ref. 6) 

externally applied load to beam AB, lb or in- lb 

linear torispherical coordinates of sphere, of torus, and of torus 
reference ring, respectively, in. 

R value at joint centroid, R = AR + R„ (fig. 6), in. 

£> 

radius, in. 

cross-sectional area, in^ 
temperature, °F 

zero stress temperature of composite structure, °F 
thickness, in. 

displacement in R, cp, and 9 directions, respectively, in. 
rotation about R, <p, and 9 axes, respectively, rad 
column vectors of internal displacements at end A and at end B 
column vectors of applied displacements at end B and at end A 
group of joint displacement components u^ to u g (ref. 6) 

coefficient of thermal expansion, in/in/°F 

mechanically equivalent thermal load parameters (defined by 
eqs. 18 and 19) 

angular torispherical coordinate, deg 
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V 


Poisson ratio 


3 

p density of core material, lb/in 

c 

3 

p’ gross core density, lb/in 

c 

o meridional stress, psi 

cp angular torispherical coordinate, deg 


Subscripts: 

a 

b 

c 

e 

fp 

i 

j 

n 

R, cp, 9 
s 


ablator 

bond 

core (moduli based on gross cross section) 

equivalent 

faceplate 

the ith laminate 

joint 

the nth laminate 
torispherical coordinates 
substrate 


THEORY 


Structural Idealization 

It was assumed that plane stress conditions existed in the class of laminated 
shells under investigation. The heat-protective material had a thermal gradient 
through its thickness. The temperature-dependent mechanical properties of the abla- 
tor were represented by multiple laminates, each at a constant temperature (fig. 2). 
Additional laminates were used to represent a bond layer and the elements of the sand- 
wich construction. The temperature gradient through the shell thickness was repre- 
sented in a step-function manner in which each of the laminates of figure 2 was at a 
specific temperature (fig. 3). The shell structure was idealized by a network of curved 
beams (fig. 4). 
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Figure 2. - Exploded view 
of composite structure. 



Figure 3. - Temperature distribution 
with step functions in the R direc- 
tion. 


The torispherical shell, idealized in 
figure 4, is the general structural shape 
used in the theoretical development of the 
present research; however, a similar ap- 
proach could be followed for other shapes. 
The theoretical development of the stiff- 
ness matrix for the typical beam element 
is presented in appendix A, and the neces- 
sary transformations appear in appendix B 


Coordinate System 

A torispherical coordinate system 
was used to define the geometry of the 
structural model shown in figure 4. In the 
spherical coordinate system, a point in 
space is described by its distance and direction from a fixed reference point. The 
coordinates are R, <p, and 9, which represent a distance and two angles, respec- 
tively. Torispherical coordinates have the same two surface coordinates <p and 9, 
but the coordinate R must be redefined as R = f(Rj, Rg), in which Rj is the torus ra- 
dius measured from the torus reference ring and Rg is the distance from the refer- 
ence point to the loci of points describing the torus reference ring. If R^ is set equal 

to zero, the torispherical coordinate system is reduced to the more familiar spherical 
coordinate system. 



Figure 4. - Structural idealization of 
the Apollo aft heat shield. 
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Procedure 


A stiffness matrix method was followed. In the formulation of the procedure, it 
was necessary to obtain the stiffness matrix of the typical structural element. Since 
a flexibility matrix was easier to obtain, the stiffness matrix was determined by the 
inversion of the appropriate flexibility matrix (ref. 6). In matrix notation, this may 

be expressed as k BB = f BB ~ The flexibility matrix was assembled through the ap- 
plication of the method of consistent distortion (see fig. 5 and appendix A for the devel- 
opment of the flexibility matrix f BB ). In general, a displacement at one end of the 


NASA-S-66-11549 DEC 6 



Figure 5. - Unit loadings for influence coefficients to type 1 beam elements. 


typical structural element will cause actions to be realized at both ends of the element. 
These actions are shear force normal to the shell midsurface, axial force tangent to 
the shell midsurface, shear force tangent to the shell midsurface, bending moment nor- 
mal to the shell midsurface, torsional moment, and bending moment tangent to the 
shell midsurface. If the typical structural element AB, fixed at end A, undergoes a 
unit displacement at end B, six actions occur at end A and six actions occur at end B. 
Thus, k BB is obtained through successive unit displacements at end B while end A 

is held in a fixed condition. A coupling stiffness matrix k^ B relates the actions at 

end A due to unit displacements at end B. The total stiffness matrix for a member 
AB is 


k AB - 


k AA k AB 

lr Ir 

L BA BB 


( 1 ) 
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AB 

The beam element stiffness matrix k is of the 12th order, 
trix A^ b was used to express the submatrices k AB’ k BA’ 

of k-- (ref. 6). Thus 

Jt>J3 


A transformation ma- 
and k. . as a function 


k AB “ A AB k BB 


k BA k AB 


k AA _ A AB k BB A AB 


( 2 ) 


The transfor ma tion matrix A^ B is developed in reference 7. The actions applied to 
the joint by the beams are shown in figure 6. The action-displacement matrix is 
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t \ 





r \ 

pa 

> = 

k AA 

k AB 

< 

U A 

> 

P B 


k BA 

k BB 


U B 

J 





S, > 


Figure 6. - Transformation matrix 
geometry. 

stiffnesses of the beams entering the joint, 
sented by u. and the externally applied loads by p 


(3) 


The stiffness matrix method is char- 
acterized by equilibrium equations ex- 
pressed in terms of unknown displacements. 
The equilibrium equations of this method 
are expressions of joint equilibrium. The 
joints are loaded by the beam elements with 
force expressions of the form of equation (3). 
External loads are introduced directly to 
the joints. Thermal loads are introduced to 
the joints in terms of mechanically equiva- 
lent loads. The thermal load manipulation 
will be discussed later in this section. In 
matrix notation, a joint equilibrium equation 
appears as km + p Q = 0; where k is the 

joint stiffness matrix, that is, the sum of the 
The displacements of the joints are repre- 


The typical joint equilibrium equation involved only orthogonal members. Each 
joint had an orientation R, cp, 9, where R was the value of the R coordinate at the 
joint centroid. A joint and all type 1 and type 3 members entering the joint had a com- 
mon coordinate system, but type 2 and type 4 beam element data had to be rotated 
through the angles 9 = 90° and <p = <Py respectively. Type 1 and type 3 beam 
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elements were oriented with joint <p coordinate direction, and type 2 and type 4 beam 
elements were oriented with the joint 9 coordinate direction. A transformation of 
type 2 beam element actions was required so that the actions might be expressed as 
functions of the joint coordinates. A similar relationship was true between type 3 and 
type 4 beam elements. Joints occurring on the outer toroidal boundary were loaded by 
three members instead of four._ The joint in which the most complex transformation 
arose was the joint located at R, 0, 0. 


The joint at R, 0, 0 was unique in that only type 1 nonorthogonal members inter- 
sected it. The joint orientation was taken as R, 0, 0. If p. represents the action 

components of a member in terms of the joint coordinate system and 9 is the member 
orientation relative to the joint, then 


P’l = Pj 

P 2 = p 2 C0S 9 " P3 sin 0 
P 3 = P 2 sin 9 + p 3 cos 9 
P 4 = P 4 

Pg = Pg cos 9 - pg sin 9 
p'g = Pg sin 9 + Pg cos 9 


( 4 ) 


The above relationships may be expressed more concisely as functions of the transfor- 
mations discussed in appendix B: 



The forces at end A of the beam elements were rotated, and the forces at end B 
were premultiplied by the identity matrix only. 

The assembly of the joint equilibrium equations was accomplished by the premul- 
tiplication of the matrix of member stiffnesses (expanded stiffness matrix) by a connec- 
tion matrix. The matrix of member stiffnesses postmultiplied by the displacement 
vector yielded the actions on the beam. The use of a connection matrix is discussed 
in reference 6. The order of the resulting square matrix was reduced through the ap- 
plication of boundary conditions. 
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Two sources existed for the reduction of the unknowns prior to the solution of 
the set of joint equilibrium equations. These sources were boundary conditions and 
symmetry conditions. The displacement functions which describe an external support 
are fully described input quantities. For example, complete fixity of a joint may be 
realized by setting Uj through u g = 0. A joint hinged along a circumferential line is 

known to have zero displacements in the u 1 through u c directions. Settlement of the 

external supports, that is, u ^ 0, is also a valid boundary condition. If elastic exter- 
nal supports are encountered, the spring rates of the supports would be required. The 
displacement u. in question would then be associated with the spring rates in the joint 

equilibrium equations. 

Symmetry conditions may be separated into the following: (1) the special case of 
the joint at R, 0, 0, and (2) the remaining joints on the diameter defining the mirror 
image boundary, that is, when 0 = 0°. The mirror image symmetry was the result of 
a geometrical orientation of the capsule relative to the flight path in a manner similar 
to the orientation of an airplane to its flight path. In the special case (1), each beam of 
orientation (R, cp, 0) occurred in company with a corresponding beam of orientation 
(R, cp, -d). The displacements at the plane of symmetry for both cases are 

0 ' 

u 2 ^° 
u„ = 0 

> ( 6 ) 

u 4 =° f 


V° ) 


Mechanical loads were introduced directly into the joint equilibrium equations. 
Distributed load systems were approximated by a set of concentrated loads applied 
at the joints. Equivalent joint loads were calculated from the reactions of die loaded 
fixed -ended beam elements. 

The thermal gradient should be considered before discussing the calculation of 
the mechanically equivalent thermal loads. The temperatures of reentry vehicles of 
the type considered have the most significant gradient in the thickness direction of the 
ablator. This gradient may be approximated by dividing the ablator into several 
layers, each at a constant temperature, so that a step function results. It had been 
previously determined that the variation of temperature through the ablator thickness 
is parabolic, while the substrate temperature is essentially constant. The temperature 
of ablation is 1000° F. 
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The temperature of ablation determines the maximum temperature during re- 
entry. This temperature is not exceeded because of the formation of a char on the 
outer surface and because of the emission of cooling gases which limit the maximum 
temperature. 

Figure 3 illustrates the temperature distribution through the thickness as repre- 
sented by the approximating step function. The temperature in the ith layer of n lay- 
ers is given by 


T. = T g (<p, Q) + 


1000 


T s ( H(s) 


( 7 ) 


The application of thermal loads was accomplished through a set of mechanically 
equivalent thermal loads. The loads were equal in magnitude but opposite in sense 
from the fixed-ended reactions of the thermally loaded beam element. These mechan- 
ically equivalent thermal loads were designated 


and 


^AlO^ 


• A20 


■AO 


-< 


• A30 


P A40 

P A50 

l P A 60 j 




P 




B20 


BO 


= 


P 


B30 


B40 




P B50 

^ P B 60 J 


( 8 ) 


( 9 ) 
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The mechanically equivalent thermal loads at a joint were the sum of the fixed- 
ended actions of all the thermally loaded curved beams entering the joint. The solution 
of a thermally loaded curved beam was divided into two load conditions; (1) that due 
to uniform thermal strain, and (2) the nonuniform strain associated with thermal mo- 
ment. The thermal moment at end B due to nonuniform thermal strain of type 1 and 
type 3 beams was designated 


»B60 = Vl 


( 10 ) 


The restraint necessary to fix the ends of a type 1 or type 3 curved beam under uni- 
form thermal strain is developed as follows: if end A is fixed, the thermal transla- 
tions at the free end B (fig. 6) are 


U B10 = 

R^l - 

U B20 = 

R sin 

U B30 = 

0 

U B40 = 

0 

U B50 = 

0 

U B60 = 

0 


C0S V Ab) V 2 
<W2 


\ 




} 


(ID 


The compatibility equations may be written through the use of the previously developed 
element flexibility submatrix fgg and the above thermal displacements. 


f BB P B0 + U B0 0 


( 12 ) 


Solving for the desired thermal restraints at end A, it is found that 


P B0 = - k BB u B0 


(13) 
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I 




The moment b^y^, described in equation (10) for the nonuniform thermal strain, must 
be superimposed on the above P B 0 q value which was necessary to maintain the condi- 
tions of fixity under a uniform thermal strain. The loads at end A of type 1 and type 3 
beams are related to those defined for end B as follows: 
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The mechanically equivalent thermal loads from type 2 and type 4 beams may be 
similarly defined in terms of the joint coordinates as 


'W 


f 1 - COS 0 AB > 


0 "I 

P B20 


Sin 0 AB 


0 

P B30 

\ R > / 2 A B , B" k BB^ 

0 


0 

> 

P B40 


0 


-cos <p^ 

P B50 


0 


0 

< P B60 j 


^ 0 > 


sin cp. 


where R is R sin cp . for type 2 and type 4 beams. The thermal loads were developed 

in the coordinate system of figure 7 and then transformed through premultiplication 

by A_,_. The mechanically equivalent actions occurring at end A of the beam are 
B B 
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Figure 7. - Joint and type 2 coor- 
dinate relationship. 


'e A io' 


f p bi(P 

P A20 


P B20 

J P A30 

> = < 

_P B30 

P A40 


_P B40 

P A50 


" P B50 

^ P A60> 


P B60> 


The definition of the mechanically equivalent thermal loads has involved the use 
of two as yet undefined parameters. They are 





2 


a E 
a a 





+ t b a b E b 



T, ^ + t, a E (2AR + t Vt -t\ 
b / fp s s V sj\ o s ) 


(18) 
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and 


y„ = > — a E (t - T.') + t E, a, (t - T. ) + 2t f E a (t - T ) 

2 Z_j n a a \ o 1/ obb\o b/ fpss\o s ) 

i=l 

fk \ + *b E b + 2t t P E s 


(19) 


It should be noted that all laminates were temperature -dependent with respect to 
their mechanical properties and the coefficient of thermal expansion. The temperature 
of the ith laminate of n laminates has been previously described in this section. With 
reference to figure 6, it may be observed that the thermal restraint forces PgQ had 

to be supplemented by a set of thermal shear restraint forces p^Q (p,^, P40’ anc * 

PgO are equal to zero) in order to satisfy equilibrium. The set of thermal restraint 

forces for any beam element formed a self -equilibrated force system. The thermal 
moments PgQ and PgQ were obtained by summing the product of the axial restraint 

force on each laminate and the distance to the neutral axis. 


EXAMPLE APPLICATIONS OF THE PROCEDURE 


Solutions to example problems solved by the method presented in this report 
were compared with other work. A limited amount of test data was obtained on a tori- 
spherical shell which had been statically tested (ref. 8). A laminated spherical cap 
under cold-soak conditions was subjected to a closed-form analytical solution (ref. 9). 

A laminated spherical cap and a laminated torispherical shell were analyzed under 
various mechanical and thermal load conditions (ref. 9). Many grid patterns (for ex- 
ample, figs. 8 and 9) were used in the present research to establish convergent solu- 
tions. The specific example problems which were used in the correlation program 
are outlined in the following paragraphs. The results are discussed in the next section. 


Example 1 - Shallow Spherical Cap 
Under Distributed Loading 

The structure for this example was a shallow spherical segment of the Apollo aft 
heat shield which had a radius of 175. 6 inches. The structure was composed of two 
0. 008 -inch stainless -steel face sheets brazed to a stainless-steel honeycomb core. 

The core was 3/16-inch square cell of 0. 0010-inch foil gage. 


A static test was conducted in which a distributed loading was simulated over a 
torispherical shell that represented the Apollo aft heat-shield substructure (ref. 8). 
The loading simulated the maximum dynamic air loading during the launch phase. It 
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Figure 8. - Apollo aft heat-shield analysis 
using a 12, 6 grid pattern. 



Figure 9. - Apollo aft heat-shield analysis 
using a 9, 4 grid pattern with the struc- 
ture fixed at the bolt circle. 


was assumed that the structure was fixed 
at the bolt circle (fig. 10). The bolt circle 
thus formed the boundary of a spherical 
cap. This spherical cap under distributed 
loading was analyzed using a finite- 
difference technique. An analysis of the 
same structure was performed under iden- 
tical conditions using a finite -element 
technique (ref. 9). The two analytical so- 
lutions and the static test results are com- 
pared with results obtained from the 
method presented in this report. 


Example 2 - Torispherical Shell 
Under Distributed Loading 

The complete torispherical structure 
used in Apollo Test Request 209 (ATR-209, 
ref. 8) was supported at the bolt circle and 
at the torus. It was the same sandwich 
structure as was used in example 1. An 
analysis (ref. 9) of this structure under the 
ATR-209 loading condition was used for 
analytic correlation. The static test re- 
sults are also included. 



Figure 10. - ATR-209 pressure 
loading. 
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Example 3 - Torispherical Shell Under Water Impact 

During the water landing of an Apollo command module, the torispherical shell 
is subjected to a pressure loading over a limited area. The extent of the area is a 
function of time (figs. 11 and 12), and is defined by the area of water in contact with the 



Figure 11. - Water impact landing, 
contact angle cp. 


Figure 12. - Water impact landing pres- 
sure. 


spacecraft at any instant of time. The problem has many dynamic ramifications af- 
fecting the loading, stability, and structural response. This investigation was limited 
to the elastic analysis of the quasi-static structure. 


Example 4 - Shallow Spherical Shell at -150° F 

The same structure was used in example 4 as in example 1, but an ablator with 
a constant thickness of 1. 84 inches was considered to be bonded to the substrate. The 
structure was subjected to a cold-soak (-150° F) thermal loading. A closed-form so- 
lution to this problem was reported in reference 9. Gallagher (ref. 9) also used his 
finite-element method to analyze the same structure under identical conditions. Both 
the closed-form and numerical-technique solutions were used for correlation purposes 
with results obtained from the present method. 


Example 5 - Torispherical Shell at -150° F and at Reentry 

The torispherical shell of example 2, with a variable thickness ablator, was 
exposed to two thermal load environments. The first was under the cold-soak (-150° F) 
conditions of example 4, which are representative of the translunar phase of spaceflight. 
The second thermal loading was a representative reentry hyperthermal load condition. 

In the performance of its function, the ablator limits the exterior temperature to the 
temperature of ablation (1000° F). The stainless-steel sandwich substrate had no 
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gradient and was always less than 600° F. The gradient through the ablator is para- 
bolic from the 1000° F surface temperature to the temperature of the substrate. 
Gallagher (ref. 9) analyzed this problem using his finite-element method. The analyt- 
ical results for the two thermal load conditions using the present method are corre- 
lated with the work of reference 9. 

The spherical cap example problems and the spherical portion of the torispherical 
shells are actual representations of the substrate of the original Apollo aft heat shield. 
The torus area is composed of an outer shell backed up by corrugated ribs or "hat" 
sections. The present method was programed to accept basic laminate geometry and 
to convert these input data to the required stiffness parameters. Each laminate had 
its own moduli and mechanical properties, which were functions of material and tem- 
perature. It was convenient to introduce a structurally equivalent stainless-steel 
sandwich substrate in the torus area. 

The temperature at which the composite laminated structure of example 4 and 
example 5 experienced zero stress was 185° F. The ablator was bonded to the sub- 
strate at a somewhat higher temperature, but a relaxation of stresses occurred during 
the curing cycle of the bond. 


DISCUSSION OF RESULTS 


The research described in this paper was applied to the five example structures 
discussed in the previous section. Comparisons of the results are discussed in this 
section. 

Example 1 - Shallow Spherical Cap Under Distributed Loading 

The displacement normal to the surface of the shallow spherical shell for the fi- 
nite difference analysis, the reference 9 analysis, and the present analysis is presented 
in figure 13. The displacement values obtained agree more closely with the ATR-209 
test values than the displacement values from either reference 8 or reference 9. 

A comparison of meridional membrane stresses is presented in figure 14. The 
agreement among the analytical results is very good. Additionally, the stress at the 
apex (c o = 0) obtained from the equation for hydrostatic loading (in ref. 10, Flugge 
points out that the stress in a shallow spherical shell under hydrostatic loading ap- 

PR 

proaches the stress in areas away from the boundary) is in remarkable agreement 

at this point with the value from the present research. The differences between the 
recorded test values and the analytical values were never resolved, and some question 
remains as to the accuracy of the test values. This circumstance, combined with the 
observation that the three independent analyses are in close agreement, indicates a 
good possibility that the higher values are more reliable. The test data and the analyt- 
ical data would aline in a very satisfactory manner if the strain gage factor of 2. 06 
were applied to the test values from reference 7. This is also true of the meridional 
outside face sheet stresses shown in figure 15. It should be noted that the shape of the 
test data curve is the same as the author’ s analysis curve and that the 2. 06 gage factor 
would again provide the desired alinement. 
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meridional outside face sheet stress, pst x 10 



Figure 13. - ATR-209, displacement 

along axis of symmetry, fixed at bolt 
circle (unpublished data refers to a 
finite difference method). 


Figure 14. - ATR-209, comparison of 
meridional membrane stresses along 
diameter of symmetry (unpublished data 
refers to a finite difference method). 


Example 2 - Torispherical Shell 
Under Distributed Loading 

This example problem served pri- 
marily to demonstrate the convergence of 
the stiffness matrix solution presented in 
this paper. Many different grid patterns 
and grid densities were investigated. The 
stresses were compared previously in ex- 
ample 1 and were not significantly differ - 
• Reference 9 ent in example 2. Thus, in reference 7 

■ Present research the convergence was discussed in terms of 
0 Re,e ' enceS displacement relative to the shell surface. 

Little discernible difference could be 
found in the various grid patterns. It was 
concluded that the method would yield sat- 
Figure 15.- ATR-209, comparison of isfactory results for this structure without 

meridional outside face sheet stresses attendant convergence problems, 
along diameter of symmetry (unpub- 
lished data refers to a finite difference 
method). 
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Example 3 - Torispherical Shell Under Water Impact 


The most severe loading of the Apollo aft heat- shield substructure occurs during 
the water landing. The critical internal loading is the R coordinate shear at the bolt 
circle. The command-module attitude just prior to impact is shown in figure 11. The 
effect of pitch angle cp on the impact loading and the relationships of pressure-time- 
wetted radii are illustrated in figure 12. A pseudostatic problem was defined by se- 
lecting a particular time after impact and noting the associated pressure-wetted radius 
shown in figure 12. A pitch angle of 15° and a wetted radius of 20 inches were run in 
the example (fig. 16). This condition occurs 0. 003 second after impact, and there is 
a pressure of 150 psi over the circular area. 


Critical values of shear occurred in the core at the bolt circle. Thus, in addition 
to a comparison of displacements normal to the surface (fig. 17), a comparison of core 
shear at the bolt circle was made (fig. 18). The displacement correlation with the 



Figure 16. - Water impact area location 
relative to bolt circle of the Apollo aft 
heat shield. 



Figure 17. - Water impact, u^ dis- 
placement along axis of symmetry 
of the Apollo aft heat shield. 



Figure 18. - Water impact, shear at the 
bolt circle of the Apollo aft heat 
shield. 
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results from reference 9 was fair. A comparison of the external reactions from ref- 
erence 9 with the external reactions from the present research at the bolt circle was 
excellent. Difficulty was encountered when the internal stresses of the structural 
model were compared with those of the actual structure. 

The support at the bolt circle of the actual structure is made up of 59 bolts, each 
spaced at 6. 2 inches with a free span of less than 3 inches. The free-span distance 
from node point to node point along the bolt circle on the structural model was 
16. 4 inches. Circumferential shear occurred on the structural model which was known 
to be at variance with the actual structure. In the calculation of core shear, this 
circumferential shear was redistributed in the same proportion as the meridional core 
shear. The applied load is introduced at the node points, and a portion of the water- 
impact load must be placed at the support node points. Thus, it is necessary to intro- 
duce this portion of the core shear into the final core shear. Each joint was in 
equilibrium and this equilibrium was maintained in the resolution of the above factors 
to obtain core shear. 


Example 4 - Shallow Spherical Shell at -150° F 

A closed-form solution under cold-soak (-150° F) conditions of a shallow spheri- 
cal cap of stainless- steel sandwich substrate covered with a 1. 84-inch-thick ablator 
was reported in reference 9. The equations for the closed-form solution were derived 
in reference 9. The comparison of the displacement normal to the shell surface with 
comparable values by the present method was excellent (fig. 19). The displacement 
values of reference 9 were somewhat larger than the values obtained in the closed-form 
solution. 


Example 5 - Torispherical Shell at -150° F and at Reentry 


displacements for the torispherical shell under cold-soak 
The values are uniformly less than those of the 


The comparison of u^ 
conditions is shown in figure 20. 

NASA-S-67-536 



Uj displacement along axis of 

symmetry of the Apollo aft heat 
shield (correlation with a closed- 
form solution). 



Figure 20. - Cold-soak condition, 

Uj displacement along axis of sym- 
metry of the Apollo aft heat shield. 
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analysis of reference 9. The correlation was considered to be satisfactory and com- 
patible with the correlation shown in figure 19. The present analysis showed excellent 
agreement with the closed-form solution in figure 19. Consequently, it should be ex- 
pected that the analysis of reference 9 would indicate a slightly more flexible struc- 
ture. 


A comparison of the cp normal stresses in the outside face of the ablator is 
shown in figure 21, while the cp normal stresses in the inner face sheet of the sub- 
strate obtained from the analyses of reference 9 and of the present research are illus- 
trated in figure 22. The EaT actions are very dependent on grid size as they represent 
the fixed-ended conditions of a curved beam element under thermal loading. In the 
limit, as the grid size becomes more dense, the thermal stress would approach the 
value Ea T. It thus appears that satisfactory displacements are possible in thermal 
problems without a comparable correlation of stresses. 



• Reference 9 

■ Present research 
9,8 grid 


s r 


• Reference 9 


Figure 21. - Cold- soak condition, ablator 
normal stress in the cp direction along 
axis of symmetry of the Apollo aft heat 
shield. 


Figure 22. - Cold-soak condition, inner 
face sheet normal stress along axis of 
symmetry of the Apollo aft heat shield. 



The displacements .experienced during 
the reentry thermal loading were an order 
of magnitude less than those which occurred 
during the cold-soak (-150° F) thermal load- 
ing. Two different grid patterns were run 
to represent the analysis of the reentry 
case. The solutions were essentially iden- 
tical (fig. 23); thus, a convergent solution 
was obtained for the present analysis. The 
larger percentage differences in this case 
between the results of the analysis of ref- 
erence 9 and the present analysis are un- 
explained. 



Figure 23. - Reentry condition, Uj dis- 
placement along axis of symmetry of 
the Apollo aft heat shield. 
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RECOMMENDATIONS FOR FURTHER RESEARCH 


There are several avenues of research which may be pursued to improve the 
procedure. These may be divided into three categories: (1) refinement of the theory, 

(2) application of the procedure to a greater variety of structures to determine the lim- 
itations of the procedure, and (3) more efficient use of the digital computer. 

The coupling of the biaxial stresses has been mentioned earlier in this section. 

To effect a tractable solution, it was necessary to assume equal biaxial stresses in the 
sizing of the axial load -carrying capability of the typical beam elements. A possible 
approach to this problem would be first to assume equal stresses in the meridional and 
circumferential directions, and then to use the stresses obtained to size axial capabil- 
ity for a second solution. This process could be repeated until the desired accuracy 
was obtained. 

It has been assumed that the change in geometry due to loading of the structure 
would have negligible effect on the solution. The influence of a generalized beam- 
column effect has been neglected. The improvement in the effectiveness of the pro- 
gram would be most noticeable in structures loaded near the buckling point. Thermal 
effects in reference 7 were limited to axial restraint and to restraint of the moment 
due to a temperature gradient through the thickness. The effect of the temperature 
variation in the shell midsurface was approximated in that each beam element was at 
a temperature constant in the <p and 0 coordinates, but variable in the R coordi- 
nate. The temperature at any location on the structure was actually an approximation 
of a three-dimensional temperature variation, since the R variation in temperature 
was represented by multiple-stepped layers. Additional research could be directed 
toward more sophisticated thermal stress representation. 

In the discussion of the results, it was noted that the thermal loads at the fixed 
ends of the beam elements are very dependent on grid size. Further research is rec- 
ommended to seek a practical grid density at which this problem would diminish to 
engineering acceptance. 

Curved beam elements were used in the present study. The displacement of the 
elastic axis due to a curved-beam effect was found to be negligible for the example 
structures (ref. 11) and was not included in this study. In structures which would re- 
quire this consideration, the theory and computer programing are straightforward. 

A modified von Karman membrane analogy was used to calculate the torsional 
capability of the beam elements in the present research. It is recommended that an 
experimental study of laminated composite beams be undertaken to improve the ap- 
proach to torsional stiffness representation. 

The exploration and extension of the program will have to be preceded by a 
learning period in which experience in the use of the computer program is attained. 

The good results experienced in the use of the present method on the Apollo aft heat- 
shield substructure were dependent on good structural idealization and realistic bound- 
ary conditions. The beam network may work satisfactorily for a wide range of shell 
structures, but this can only be determined through the application of the computer 
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program to many more structures than demonstrated by the example problems. Once 
the limitations of the computer program are understood and confidence in its use has 
been established, it would be possible to extend the program into the areas of linear 
elastic stability (ref. 12) and linearized plasticity (ref. 13). Both of these areas were 
outside the scope of the study presented in this report, but they are noteworthy direc- 
tions in which further research would prove beneficial. 


CONCLUSIONS 


The stiffness matrix solution yielded highly satisfactory results for the Apollo 
aft heat-shield substructure. Success was achieved primarily because of the adequacy 
of the structural idealization. The frame networks of beam elements used to represent 
the shells of the example problems were demonstrated to be adequate for engineering 
purposes. Experience with the program will be required to establish the full range of 
structures which may be validly analyzed. 

To effect tractable solutions of complex structures, it is necessary to make en- 
gineering compromises in all methods. The effect of Poisson's ratio was not com- 
pletely accurate in the present report, but it was sufficiently accurate to obtain good 
engineering results. A plane stress condition was assumed in order that beam element 
stiffness matrices might be calculated; that is, the meridional and circumferential 
stresses were assumed equal so that the beam elements might be sized to reflect the 
Poisson effect. The use of concentrated loads to simulate a distributed load system 
could lead to difficulty if the grid size were too coarse. The assumption that the 
deformed geometry does not affect the solution ignores the effect of generalized beam- 
column action. Incorporation of this effect in linear elastic stability analysis has been 
accomplished by R. H. Gallagher and Joseph Padlog. 

The attainment of a convergent solution for the Apollo aft heat shield was easily 
accomplished with coarse grid patterns. It is not expected that an indefinite refine- 
ment of the grid would lead to the "exact" solution. The adequacy of the beam ele- 
ments to represent the idealized structure will break down in the limit. The use of 
extreme aspect ratios in the selection of grid patterns led to some difficulty in the sat- 
isfaction of equilibrium equations. This difficulty was not serious in the example prob- 
lems but should be recognized. The limitations discussed did not adversely affect the 
structural analysis of the Apollo command module aft heat shield. It is concluded that 
composite shells of similar geometry and loadings may be successfully analyzed with 
engineering accuracy. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, January 20, 1967 
914-5020-01-72 
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APPENDIX A 


FLEXIBILITY MATRIX 


The details of the development of the flexibility matrix f^^ were not essential 

to the general development in the section on theory. Therefore, appendix A was des- 
ignated for the derivation of f DD . 

DD 

Four structural elements were used in the analysis. It is emphasized that all 
elements were of the same structural type, that is, laminated beams. Each beam ele- 
ment was assumed to carry six generalized actions (three forces and three moments). 
The division into four types was necessary because of the variation in geometry of the 
example structures. The redefinition of geometry made it possible to use a single set 
of equations for all flexibility calculations. Type 1 and type 2 beams were on the 
sphere in the meridional and circumferential directions, respectively. Type 3 and 
type 4 beams were on the torus in the meridional and circumferential directions, 
respectively. 

In the determination of the structural model, it was assumed that each beam ele- 
ment represented one -half of the area on either side of the beam. It was also assumed 
that the Poisson effect could be introduced into the present method through the modifi- 
cation of the C. terms. The bending stiffness was increased by — and the axial 

1 1 " v 
stiffness by -j — -. 

The evaluation of the stiffness data is to a great extent self-explanatory; however, 
an explanation of the effect of temperature on the stiffness of the elements is in order. 
The specific treatment of the effect of temperature on the stiffness of the beam ele- 
ments in the Apollo command module aft heat-shield substructure is discussed in this 
appendix; however, the discussion is applicable to a general class of composite shell 
structures. The composite nature of the structure is shown in the exploded view of 
figure 2. The steel sandwich substructure was at a constant temperature, but the abla- 
tor temperature varied in a parabolic manner from the temperature of ablation 
(1000° F) at the outside surface to a temperature T^ at the bond. The step-function 

approximation of temperature is shown in figure 3. In the range of temperatures en- 
countered, the moduli of the steel substrate were not significantly affected. However, 
the small effect (approximately 10 -percent reduction at 600° F from room temperature 
values) was incorporated into the computer program. The mechanical properties of 
the bond and ablative layers were very temperature -dependent. Their moduli and the 
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steel substrate moduli were represented in the computer program as polynomials in 
the temperature. The equation 1 for the modulus of rigidity of the core G c was an em- 
pirical formulation, and is based on the gross cross section of the core. Hoff (ref. 14) 
presents a formal derivation of G but the results did not justify the increased com- 

plexity of its use. 

hi the calculation of the member stiffness, the composite nature of the structure 
made it necessary to associate each laminate with its correct moduli. The moduli 
were temperature-dependent and had to be considered in the calculation of the centroi- 
dal coordinate R. The R coordinate to the outside face sheet was Rg, and the cen- 

troidal R coordinate (R) was obtained by calculating the movement AR from Rg. 

The moduli and geometry for the beam elements were obtained by averaging the 
terminal values, that is, the values at the joints. A discussion of the method used to 
obtain the torsional stiffness of the typical structural element appears at the end of 
this appendix. 


Formation of Beam Element Flexibility Matrix, f 


BB 


BB 


Bl, 

Bl 

f Bl, 

B2 

f Bl, 

B3 

f Bl, 

B4 

f Bl, 

B5 

f Bl, 

B6 

B2, 

Bl 

f B2, 

B2 

f B2, 

B3 

f B2, 

B4 

f B2, 

B5 

f B2, 

B6 

B3, 

Bl 

f B3, 

B2 

f B3, 

B3 

f B3, 

B4 

f B3, 

B5 

f B3, 

B6 

B4, 

Bl 

f B4, 

B2 

f B4, 

B3 

f B4, 

B4 

f B4, 

B5 

f B4, 

B6 

B5, 

Bl 

f B5, 

B2 

f B5, 

B3 

f B5, 

B4 

f B5, 

B5 

f B5, 

B6 

B6, 

Bl 

f B6, 

B2 

f B6, 

B3 

f B6, 

B4 

f B6, 

B5 

f B6, 

B6_ 


(Al) 


The equation for the modulus of rigidity of the core was developed by North 
American Aviation. 
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It was stated in the section on theory that the flexibility matrix fgg would yield 

the stiffness matrix kgg when inverted. The elements of fgg were obtained through 

the use of the auxiliary load method (virtual work). The load diagrams resulting from 
the application of unit loads on a type 1 element AB, free at end B and fixed at end A, 
are shown in figure 5. These load expressions were inserted in the displacement equa- 
tions of reference 6. As quoted from reference 6, the displacement is as follows: 



where 

A x , A^, A z are the cross-sectional areas associated with the x, y, z actions 
dl is the incremental length along the beam 
E is the modulus of elasticity 
G is the modulus of rigidity 

I x , 1^, I z are the moments of inertia associated with the x, y, z actions 
M X is the twisting moment about the x-axis 
is the bending moment about the y-axis 

Z 

M is the bending moment about the z-axis 

N is the axial force in the x direction 

is the shear force in the y direction 
z 

N is the shear force in the z direction 
n is the unit quantity of N X , N^,N Z , M X , M^,M Z at r 
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I ii nil nn ■■ ■ n ii ■■■■■■■■ in mu 


U 


is the displacement at r 
signifies integration around the 


whole frame 


A summary of the flexibility equations appears in table A-I. The definition of the C. 
expressions appearing in the table are shown later in this appendix. 

Type 1 and Type 3 Members, Cj 


c - 1 - v 
1 “ 2EI e 

(A3 a) 

c = 1 - 1,2 
2 2ei r 

(A3b) 

r - 1 

'■'3 " 2GJ 

e 

(A3c) 

C = 1 

4 2GS r 

(A3d) 

c = 1 

^5 2GS 0 

(A3e) 

c = 1 - v 
^6 2ES 

(A3f) 

C 7 = ^AB 

(A3g) 

C 8 = Sin ^AB 

(A3h) 

C 9 = C0S ^AB 

(A3i) 

sin <p AB cos cp AB = C ? - C g C 9 

(A3j) 
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C 11 = ^AB + Sin ^AB C0S ^AB = C 7 + C 8 C 9 

(A3k) 

C 12 = 2 cos ^AB " 2 + sin2(p AB = 2c 9 - 2 + C 8 2 

(A31) 

C 13 = C 11 + 2C 7 - 4C 8 

(A3m) 

C 14 = C 11 - 2C 8 

(A3n) 

C 15 =1 - C 9 

(A3o) 

C 16 = C 8 - C 7 

(A3p) 


Type 1 and Type 3 Member Stiffness Data (See fig. 4) 

R = AR + Rg (A4) 
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b “ 
F T _ <P 

ei r‘Tt 


n t 

L, E a( T l)¥ + % + 2 Vfp 

i=l 


G = 


E 


2(1 + v) 


p'\ 1.54 

°c- 2 - 4 *lp?i 0 


GJ (See section on calculation of torsional stiffness) 
e 



i=l 


t G a( T i)T + % + 2 G s t fp 
_ i=1 

E E a( T i)^ +E bV^'fp 

1=1 -I 

Type 2 and Type 4 Members, C 

l-, 2 

C 1 = 2EI R 

1 2 

c - 1 ~ v 

2 2EI 

<P 

C = _JL_ 

°3 2GJ 

e 


GS e = b <? 


ES = b 


<p 


(A7) 

(A8) 

(A9) 

(A 10) 

(All) 
(A 12) 

(A 13a) 
(A 13b) 
(A 13c) 
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(A 13d) 


C A = 


4 ” 2GS 


9 


C c = 


5 2GS 


C* = 


R 
1 - v 


6 2ES 


C 7 = 0 AB 


C 8 = Sln 0 AB 


C 9 = C0S 0 AB 


C 10 " e AB " Sin 0 AB C0S ^AB C 7 - C 8 C 9 


C 11 = 0 AB + Sin 0 AB C0S 0 AB = C 7 + C 8 C 9 


C l2 = 2 cos » AB - 2 + sin 2 0 AB = 2C g - 2 + Cg 2 


C 13 = C 11 + 2C 7- 4C 8 


C 14 " C 11 “ 2C 8 


C 15 " 1 ' C 9 


C 16 “ C 8 “ C 7 


(A13e) 

(A13f) 

(A13g) 

(A13h) 

(A13i) 

(A13j) 

(A 13k) 

(A 131) 

(A13m) 

(A13n) 

(A13o) 

(A13p) 


31 



Type 2 and Type 4 Member Stiffness Data (See fig. 4) 


EI R = TT| 


n 

V 


+E b t b + 2E s t f P 
i=l 









+ AR - 



Note: Use AR as defined for type 1 beams. 


G = 


2(1 + v) 


GJ e (See section on the calculation of torsional stiffness) 


GS 


<P 



G b*b 


2G t. 
s fp 


GS„ 

R 





2G t, 
s fp 


ES = b 


n t 


(A14) 


(A15) 


(A16) 


(A17) 


(A18) 


(A19) 
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Calculation of Torsional Stiffness 

The conventional torsion formula is limited to circular cross sections. Prandtl 
developed a membrane analogy theory (ref. 15) for noncircular cross sections. This 
enables a valid use of the conventional torsion formula with the introduction of a geo- 
metric factor J g . The modulus of rigidity G is constant throughout the cross section 

in the Prandtl theory. In the present application, the typical structural element was 
composed of a laminated rectangular cross section, with each laminate possessing an 
arbitrary modulus of rigidity. To account for the different moduli of the laminates, 
the modulus of rigidity of the gross cross section was taken to be proportional to (1) the 
laminate cross-sectional area, (2) the laminate modulus of rigidity, and (3) the distance 
from the center of twist in the direction of the shortest dimension of the cross section. 
This was compatible with the assumption of an equal angle of twist for all laminates in 
a cross section, and of a linear shear stress distribution across the shortest dimension 
(ref. 16). The foregoing may be expressed mathematically as 


G 

e 



G.t.d. 

ill 



(A20) 


and its application to the Apollo command-module aft heat shield is found in figure A-l. 

The use of the equivalent modulus of rigidity 
G g in the modified Prandtl theory follows. 

In reference 15 


M t = k 1 G0(2a) 3 2b (A21) 


or in the notation of the present research 


L, d„ G., Id7, 

Laminate i ■ 1 ' 1 • 

in. in. psi in 2 

Ablator 3 1 1 1.728 100 000 1.728 

Bond 2 0.020 0.718 414 000 0.01436 

Faceplate 3 0.008 0.704 12 000 000 0.005632 

Core 4 1.984 0.292 27 500 0.579328 

Faceplate 5 0,008 1.288 12 000 000 0.010304 

Z — * - - - 2.337624 

a Ablator may be multilayered For temperature dependence. 


Figure A-l. - Example calculation 
of torsional stiffness. 
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I 


I III II Ill I 


A comparison with the well-known torsional formula cp = n T reveals G J 

uJ e e 

e 

to equal ^ t.. The constant is taken from reference 15 and is tabulated 


below. 


b y t. 
1 

k l 

1.0 

0. 1406 

1.2 

0.166 

1.5 

0. 196 

2.0 

0.229 

2.5 

0.249 

3 

0.263 

4 

0.281 

5 

0.291 

10 

0.312 

OO 

0.333 
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TABLE A- 1. 

- SUMMARY OF DISPLACEMENT INFLUENCE COEFFICIENTS 


[r used for type 1 and type 3 beams I 

j_R sin <p used for type 2 and type 4 beam^ 


1 

rc 6 c io 

2 

RC 6 c 8 2 

3 

4 

5 

6 

+RC^C^ j 

+(R) 3 c 1 c 10 

RC 6 C 8 2 
— 2 

-rc 4 c 8 2 

+(r) 3 c 1 c 12 

rc 6 c h 

0 

0 

0 

-2(R) 2 C 1 C 

9 

- RC 4 C 8 

+ (R) 3 CiCi2 

+rC 4 C 10 

♦® 3 c,c 13 

0 

0 

0 

-2(R) CjC 



2 rc 5 c 7 

-(R) 2 C 3 Ci2 

-(R) 2 C 3 C 14 


0 

0 

+ (R) 3 C 3 C 13 

+ (r) 3 c 2 c io 

+ (r) 2 c 2 c 8 2 

-(R) 2 C 2 Cio 

0 

0 

0 

-® 2 c 3 c 12 

*®) 2 °2 C 8 2 

RC 3 C 10 

+RC 2 C 11 

SC 3 C 8 2 

- SC 2 C 8 2 

0 




— 2 



0 

0 

-(R) 2 C 3 c i4 

RC 3 C 8 

SC 3 C H 

0 



-(r) 2 c 2 Cio 

' 

- RC 2 C 8 2 

+RC 2 C 10 


■2(R) 2 C, C, _ 

“2(R) 2 C,C, « 

0 

0 

0 

2RC^C n 





APPENDIX B 


TRANSFORMATION MATRICES 


A congruent transformation was used to obtain the total member stiffness matrix 
once the stiffness submatrix k-,,-, was obtained from the inversion of the flexibility 

matrix A typical type 1 beam element is shown in figure 6. A system of dis- 

J3J3 

placements u g induced actions p BB at end B and p^ B at end A, such that 


P BB k BB u B 


(Bl) 


P AB " k AB u B 


(B2) 


The actions p BB were expressed in the axes at end B, while actions p^ B were ex- 
pressed in the axes at end A. These two sets of actions formed a system of equilib- 
rium equations acting on the member AB. If the actions p BB are transformed to the 

axes at end A, the equilibrium equations are 


p AB + A AB p BB = 0 


(B3) 


and thus 


P AB " A AB P BB 


(B4) 


Therefore 


k AB U B " A AB k BB U B 


(B5) 


and 


k AB " A AB k BB 


(B6) 
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Similarly 


k BA " _A BA k AA 


(B7) 


The stiffnesses and k^g are also related to one another, 

ciprocal law it is known that = k^g, thus 


From the re- 


" A BA k AA “ 



(B8) 


Since the transpose of the product of two matrices is equal to the product of their 
transposes in reverse order 


and 


Since 


and 


_A BA k AA “ " k BB A AB 


k AA “ A BA k BB A AB 


A - A 

A BA " a AB 


k - k’ 
K BB *BB 


k AA " A AB k BB A AB 


(B9) 

(BIO) 

(Bll) 

(B12) 

(B13) 


The above submatrix stiffness definitions may be collected into a single matrix 
expression for the actions on the member AB in terms of the end displacements u^ 

and Ug. 
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PA ) | k AA k AB| f u A 


P B k BA k BB U B 


(B14a) 


P a) A AB k BB A AB _A AB k BB f U A 


P B I “ k BB A AB 


(B14b) 


A similar approach was used in the calculation of f BB for type 2 beam elements. 

An additional congruent transformation was required to orient the forces and displace- 
ments with the joint orientation. The relationship of the coordinates in the type 2 beam 
formulas to the joint coordinates is shown in figure 7. The transformation matrix 
manipulation follows. 


p A k AA k AB |A 


P B k BA k BB U B 


(B15) 


A A'A 0 k AA k AB ( U A 


0 A B'B k BA k BB U B 


(B16) 


P A' ) A A'A k AA A A'A k AB f U A 


A B'B k BA A B'B k BB l U B 


(B17) 


A A'A k AA A A'A k AB | A A’A 0 


A B’B k BA A B’B k BB‘ 1° A DtTDM u n* 


B’B B 


(B18) 


l A'A k AA A A'A A A'A k AB A B*B U A’ 


B'B BA A’ A B’B BB B’B I B* 


(B19) 
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The type 1 beams which meet at the apex required an additional congruent trans- 
formation of actions and displacements to make each compatible with the single set of 
joint coordinates at zero. The zero joint is oriented such that the local cp coordinate 
makes an angle of 0° with the 6 coordinate of the fixed reference point of the shell. 
This is illustrated as 


Pi 


B 


k AA k AB 


k BA k BB 


|u. 


u T 


(B20) 


Transform actions at end A to the zero joint orientation 


0A 


B 


A 0A 0 


k AA k AB 


k BA k BB 


|u, 


U T 


(B21) 


'oa! 

A 0A k AA 

A 0A k AB 



k BA. 

k BB_ 

W 


(B22) 


Then transform the displacements at end A to the zero joint orientation 
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The calculation of the thermal load vector required rotational transformations 
for all type 2 beam elements and those type 1 beam elements common to the apex. 
The use of the required transformations is as follows: 


Type 1 beam elements 



(B25) 


Type 2 beam elements 



(B26) 


The transformations discussed above required the use of special transformation 
matrices. These matrices are developed in detail in reference 7. 
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